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ABSTRACT 

We report the discovery by the Gemini NICI Planet-Finding Campaign of two low- 
mass companions to the young AOV star HD 1160 at projected separations of 81 ± 5 
AU (HD 1160 B) and 533 ± 25 AU (HD 1160 C). VLT images of the system taken 
over a decade for the purpose of using HD 1160 A as a photometric calibrator confirm 
that both companions are physically associated. By comparing the system to members 
of young moving groups and open clusters with well-established ages, we estimate an 
age of 50/40 HD 1160 ABC. While the UVW motion of the system does not 

match any known moving group, the small magnitude of the space velocity is consistent 
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with youth. Near-IR spectroscopy shows HD 1160 C to be an M3.5 it 0.5 star with 
an estimated mass of 0.22lQg4 Mq, while NIR photometry of HD 1160 B suggests a 
brown dwarf with a mass of Mj^p. The very small mass ratio (0.014) between 

the A and B components of the system is rare for A star binaries, and would represent 
a planetary-mass companion were HD 1160 A to be slightly less massive than the Sun. 

Subject headings: brown dwarfs — instrumentation: adaptive optics — planetary sys- 
tems — planets and satellites: detection — stars: individual (HD 1160) 



Introduction 



There are almost 700 extrasolar pl anets discovered to date, most of thein at small separation s 



(<5 AU) from their parent stars (e.g. iButler et al. 



2006 



Marcy et al. 



2008 



Mavor et al.l l2009l ). 



In addition, 1235 tr ansiting planet candidates have been discovered by NASA's Kepler Mission 
(jBorucki et al.ll201ll ). effectively tripling the census of exoplanet candidates. Unlike radial velocity 
and transit surveys, which are responsible for the vast majority of these discoveries, direct imaging 
allows us to access planets at much larger separations (> 5 A U) . Direct imaging surveys for gas-giant 



planets around young stars initially retur ned null results (jMasciadri et al. 



2005 



Lafreniere et al 



20071 : iBiller et al.l 120071 : iNielsen et al.ll2008l ). and a detailed analysis of these null r esults shows that 
giant planets are rare at large (> 65 AU) separations around solar-type stars (jNielsen &: Close 
2010l ). Despite the rarity of the se objects, discoveries of directly imaged planets at large separa- 



tions have now been announced (Maro 



Lagrange et al.ll2009l : iKraus Sz Ireland 



of this new population of planets (e.g. 



set 



al.l l2008l . l20inl : iKalas et al.ll2008l : iLafreniere et al.ll2008l : 



2OI2I'). allowing the study of the comp osition and formatio n 



Bowler et al.ll20ld : ICurrie et al.ll201ll : iBarman et al.ll201ll ). 



High-mass stars ha ve recently become ex citing targets for direct imaging surve ys. Both brown 



dwarfs (e.g. HR 7329 B [Lowrance et al.ll2000l | and HR 6037 B 



Hue 



mass companions (HR 879 9 bcde, Fomalhaut b, and Beta Pic b 



amo et al. 


201c n and olanetarv 


Marois et al. 


2010: 


Kalas et al. 



20081 : iLagrange et al.l 1200911') ha ve been detected around A-type stars. Analysis of radial-velocity 



results by I Johnson et al.l (I2OIOI ) has shown that short-period (a < 2.5 AU) giant planets are more 
frequent around high-mass A stars than low-mass M stars, and the direct detection of planets around 
multiple A stars (where the higher contrast ratio should make planet de tection by direct imagin g 



more difficult) suggests this trend may continue to larger separations (e.g. ICrepp &: Johnsonll201ll ) 



The Gemini NICI Planet-Finding Campaign is a 3-year survey begun in December 2008 t arget- 
ing about 300 young, nearby stars to directly image extrasolar giant planets (jLiu et al.ll2010bl ) . The 
goals of the Campaign are to detect new planets and to investigate the properties and formation 
mechanisi n of planets at >10 AU separations. Design ed specifically to detect planets, the NICI in- 
strument (iToomey &: Ftaclasll2003l : IChun et al.ll2008l ) consists of an 85-element curvature adaptive 
optics system, a Lyot coronagraph, and t wo cameras that can image simultaneously for spectral dif- 
ference imaging (SDI: lRacine et al.lll999l ). observing through two moderate-band filters around the 
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H-hand methane feature seen in cool substellar objects ([Burrows et al.l 119971 : iBaraffe et al.ll2003l ). 
NICI also uses angular differential imaging (ADI), w hich holds the telescope rotat or fixed, allowing 
speckle noise to be separated from real companions (|Liull2004l : iMarois et al.ll2006l ). 



PZ Tel (jBiller et alJl2nim and CD-35 2722 (jWahhai et alJ 12011 



We previously annou nced the Campaig n discovery of two b rown dwarfs orbiting young stars 

. We present here the discovery 
of two low- mass companions to the young A star HD 1160, at ~80 AU (B) and ~530_^AU_(C) 
separa tions. HD 1160 (see Tabled]) is a well known photometric standard star from lElias et al 



(|l982l ) and was selected for the Cam paign due to its underluminous position on t he HR diagrarn , 
a sign of youth for field A stars (e.g. lJura et al.lll998l ). HD 1160 was observed by ISu et al.l (j2006l ) 



with Spitzer at 24 fim, and no excess was detected. We derive an age range for the system of SOI^q 



Myr and estimate masses for the two companions of SSl^g^ Mj„p 
(HD 1160 C). 



(HD 1160 B) and 0.22+°;°^ Mq 



2. Observations 

2.1. Gemini-South NICI: Near-Infrared Imaging 

NICI, the Near Infrared Coronagraphic Imager, is a dual-channel high contrast imager at 
Gemini South, designed specifically to directly image giant planets. The detectors in the two 
channels are 1024x1024 Aladdin II InSb, with a square 18.43" field of view and 17.94 and 17.96 
mas pixels in the two cameras. 

HD 1160 was first observe d on 2010 Octobe r 30 UT in the standard NICI Campaign observing 



modes (previously described in lBiller et al.ll2010l and lWahhaj et al.ll201ll ). In the combined angular 
and spectral difference imaging (ASDI) mode, where we are more sensitive to methane-bearing 
planets, images are taken simultaneously in the 4% bandpass methane on (1.652 ^um) and off (1.578 
fim) filters using NICI's two cameras. In the ADI mode, where all of the light is sent to one camera, 
we are more sensitive to faint companions at separations ^2". We obtained 45 1-minute exposure 
images in the ASDI mode and 20 1-minute exposure images in the ADI mode. In the ASDI mode, 
the primary is not saturated as it is observed through a partially transmissive coronagraphic mask 
at NICI's first focal plane. 

We easily detected two companions at separations of 0.78ib0.03" and 5.15ib0.03". After using 
archival VLT data to confirm that both companions are common proper motion (Sections 12.31 and 
13. ip . we observed HD 1160 again on UT November 22, 2010 in the J, H and Ks bands (filters 
are in the MKO system) in order to measure colors for the companions. Ten images were taken 
simultaneously in the H and Ks bands with exposure times of 12.16 seconds and 5 coadds per 
image. In the J-band, we obtained twenty images in single channel mode with the same exposure 
times and coadds per image. Both companions were well detected even in individual images, but 
we chose only the best 7-10 images where the primary was centered most accurately on the mask. 
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We obtained an additional NICI dataset on UT October 21, 2011 (in the narrowband methane 
filters), which provided further astrometric confirmation. 



2.2. Keck-II NIRC2: l! and Ms band imaging 



We observed HD 1160 in the l! (3.426-4.126 //m) and Ms (4.549-4.790 ^m) bands with the 
Keck II Telescope adaptive optics system and facility near-IR camera NIRC2 on UT November 27, 
2010 without any coronagraphic mask. We used the narrow camera which has 9.961 mas pixels 
(jGhez et al.ll2008l ). In the direct imaging mode (non-ADI), we obtained 37 images at L' band, each 
with an exposure time of 0.15 seconds and 200 coadds. In the Ms band, we obtained 29 images 
with an exposure times of 0.2 seconds and 50 coadds. For each band, we used 5 dither positions 
with separations of 1". The median combination of the unshifted images was used to make sky 
images for subtraction. Both companions were clearly detected and the primary was unsaturated 
in both bands . In ad dition to the plate scale, we also use the Keck NIRC2 field rotation derived 
bv lGhez et al.l (|2008l ) to determine the astrometry of the two companions. 



2.3. VLT NACO and ISAAC: Ks, L, and L' imaging 

We examined archival data of HD 1160 taken at the VLT between 2002 and 2012 in order to 
assess whether HD 1160 B and C are common proper motion companions to HD 1160 A. These 
datasets include Ks and V (3.49-4.11 /im) adaptive optics images taken with the NACO system, 
and ISAAC Ks and L (3.49-4.07 fim) images, originally obtained for photometric calibration of 
observations of other science targets. NACO images taken in the Ks-hand have 27.06 mas pixels, 
while V images have 27.12 mas pixels. ISAAC has a Ks pixel size of 148 mas and L-band images 
have a pixel size of 71 mas. Despite their relatively short integration times (of order one minute 
per epoch), HD 1160 C (~5" separation) is cleanly detected in all the data. HD 1160 B (~0.8" 
separation) is only detected in the adaptive optics datasets and is not resolved in the seeing-limited 
ISAAC data. Position angles measured with the NACO inst rument are corrected to account for a 



-0.6 degree field rotation measured by lBergfors et al.l ()201 



Figure [T] shows the detection of both companions at each VLT NACO, Gemini NICI, and Keck 
NIRC2 epoch. Figure [2] shows the detections of only HD 1160 C at each of the seven ISAAC epochs. 
Astrometry for all epochs is presented in Table [2j 



2.4. Magellan MIKE: Optical Spectroscopy 



We acquired a spectrum of HD 1160 A and C on UT 30 December 2010 with the Magellan 
Inamori Kyocera Echelle (MIKE) spectrograph at the Magellan Clay telescope at Las Campanas 
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Observatory in Chile. We used the 0.5" sHt whi ch produces a spectral resolution of ~35,000. These 
data were reduced using the facility pipeline (jKelsonl |2003| ) . We do not attempt to recover the 
spectrum of the much closer companion, HD 1160 B, with these observations. 

Each stellar exposure is bias-subtracted and flat-fielded for pixel-to-pixel sensitivity variations. 
After optimal extraction, the 1-D spectra are wavelength calibrated with a thorium-argon arc. To 
correct for instrumental drifts, we used the telluric molecular oxygen A band (from 7620-7660 A) 
which aligns the MIKE spectra to 40 m s~^, after which we corrected for the heliocentric velocity. 
The final spectra are of moderate S/N reaching ~ 50 per pixel at 8700 A. We also observed GJ 908 
(Spectral Type Ml) as a radial velocity (RV) standard. 

To measure the RV of HD 1160 C, we cross-correlated each of 9 orders between 7000 and 9000 A 
(excludi ng those with stro ng telluric absorption) with the spectrum of GJ 908 using IRAF'f0 fxcor 
routine (jFitzpatricklll993). We use GJ 908's -71.147 km s"^ RV published by the California and 
Carnegie Planet Search (jNidever et al.ll2002l ). The zero-point of the absolute RV is uncertain at 
the 0.4 km s~^ level. We measured the RV of HD 1160 C from the gaussian peak fitted to the 
cross-correlation function (CCF) of each order and adopt the average RV of all orders with the 
standard deviation of the individual measurements. We measure the radial velocity of HD 1160 C 
to be 12.6 ± 0.3 km/s. We also observed HD 1160 A itself with the MIKE spectrograph. Given 
the few spectral lines available we were unable to get as precise an RV measurement as we did for 
the M dwarf HD 1160 C. However, measuring the RV of several lines in HD 1160 A and taking the 
average, we found that the RV of HD 1160 A is 14 ± 2 km/s, consistent with that of HD 1160 C. 



2.5. IRTF SpeX: Near-IR Spectroscopy 

We obtained moderate-resolution (R w 2000) near-IR spectra of HD 1160 A and C on 29 June 2011 UT 
from NASA's Infrared Telescope Facility (IRTF), located on Mauna Kea, Hawaii. Conditions were 
photometric with good seeing conditions (~0.5-0.6"), and the C companion was well-resolved from 
the primary star. We did not attempt to recover the spectrum of the much closer companion, 
HD 1160 B, with these observations. 



We used the facility near-IR spectrograph SpeX (jRavner et al.lll998l ) with the 0.3" wide slit 
in cross-dispersed mode, obtaining spectra from 0.8-2.5 /um. The slit was oriented at a position 
angle of 50° on the sky, compared to the parallactic angle of ~80° at the time, in order to minimize 
the light from the primary that entered the slit. In the pairwise-subtracted images, in the H and 
K bands, no evidence was seen of contaminating light from the primary compared to the very 
well-detected companion. Some contamination was seen in the J band. Since the target was nearly 



^IRAF (Image Reduction and Analysis Facility) is distributed by the National Optical Astronomy Observatories, 
which is operated by the Association of Universities for Research in Astronomy, Inc. (AURA) under cooperative 
agreement with the National Science Foundation. 
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overhead at the time of the observation (airmass ranging from 1.07 to 1.10), the effect of atmospheric 
dispersion on the broadband spectrum is minimal. The slit was not oriented perpendicular to the 
primary-companion axes, since this would cause the companion's spatial position in the slit to 
coincide with a local maximum in the light from the primary's halo. 

HD 1160 C was nodded along the slit in an ABBA pattern for a total on-source integration of 
16 min, with individual exposure times of 120 sec. We observed the primary star HD 1160 itself for 
telluric calibration, since it has a spectral type of AOV and was well-resolved fr om the companion . 
All spectra were re duced using version 3.4 of the SpeXtool software package (jVacca et al.l 120031 : 
Gushing et al.ll2004l ). The median S/N per pixel in the final reduced spectrum is 25. 



3. Results 

3.1. Astrometric Confirmation of HD 1160 B and C 

Our full set of astrometry from VLT, Gemini-South, and Keck is reported in Table El Figures 
E] and m show the relative positions of HD 1160 B and G with respect to A over nearly a decade 
(2002-2011), along with the expected relative motion of a distant background object, given the 
known proper motion and parallax of HD 1160 A. Errors on this astrometric track are computed 
through a Monte Garlo approach, using errors from the reference (NIGI, epoch 2010.8301) positio n, 



and errors in the proper motion and parallax of HD 1160 A from Hipparcos ([van LeeuwenI 120071 ) 



Figures [5] and [6] show the astrometry of the companions on the sky. The astrometric measure- 
ments cluster on top of each other, as expected for common proper motion companions (GPM), 
and do not follow the background track. 

To robustly determine whether each companion is more consistent with background or common 
proper motion, we compute the statistic using the equations 



and 



XCPM = 2^ -2 + -2 1 (2) 



^p,i ^PA,i 



where pobs,i and PAobs,i are the observed separation and position angle at each epoch (i), p^g^i 
and PAbg^i are the expected separation and position angle at each epoch from the background 
(bg) track, po and PAq are constant separations and position angles (corresponding to common 
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proper motion), and Cp^j and (JpA,i are the measurement uncertainties in the observed separation 
and position angle. 

HD 1160 B is easily confirmed as a co- moving companion, differing from the expected motion 
of a background object by 15 degrees in position angle. The fit to the background track across all 
epochs has reduced chi square x^=26.2 (with 22 degrees of freedom). These data are not consistent 
with the background hypothesis, which is ruled out at the P~ 0% level. Meanwhile, the fit to 
constant separation and position angle (i.e. common proper motion) produces x^=0.38 (dof=22, 
P=99.6%). The deviation of HD 1160 C's motion from the background track is less than that for 
B. Nevertheless the offsets are significantly more consistent with common proper motion (x^=0.62, 
dof=36, P=96.5%) than background (x^=4.95, dof=36, Pss 0%). Though examining Figure [6] 
shows that while there is a spread to the 19 astrometric epochs and they do not clearly overlap 
within error bars, the direction of the spread is orthogonal to the expected motion of a background 
object (and probably due to orbital motion, as we discuss below. In addition, we note that the 
radial velocity measurements of HD 1160 A and HD 1160 C are consistent within measurement 
uncertainties (see Section [2. 4p . as we would expect from physically associated companions. 

Figure H] suggests a small change in position angle over time for HD 1160 C, raising the 
possibility that we are observing orbital motion. Fitting a straight line with non-zero slope to the 
astrometry for HD 1160 C produces x^=0.42 (dof=34, P=99.8%), similar to the common proper 
motion fit, indicating that with our astrometric errors, we are unabl e to detect a dev iation from 



non-zero relative motion. Assuming a mass for an AO star of 2.2 Mq (ISiess et al.ll200d ). a mass for 
HD 1160 C of 0.2 Mq (see Section [3.6p . and a face-on circular orbit with semi-major axis of the 
projected separation of 531 AU, we would expect orbital motion of -^0.05 °/year, while the best-fit 
line to the PA motion gives —0.04 zb 0.03 °/year. So while we do not significantly detect orbital 
motion, the magnitude of such motion would be consistent with that expected for HD 1160 C. 



3.2. Spectral Typing of HD 1160 C 



Figure [7] shows the H-K and K-L' colors for HD 1160 B and C compared to field objects of M 
and L spectral types. The locations of the HD 1160 companions in this color-color diagram suggest 
spectral types of mid-M for HD 1160 C, and late M/early L for HD 1160 B. 

Figure [8] shows HD 1160 B and C on an MKO NIR color-color diagram that includes field 
stars, field L dwarfs, selected individual low-mass objects, a nd members of the Taurus star-forming 
region. We include in this plot photometry of field stars (IRavner et al.ll2009l) and fi eld L-dwarfs 



( Leggett et al. 



Rebolo et al 



Lowrance et al. 



.20101 '). as well a s the low mass ob jects AB Pic b (IChauvin et al.ll2005l ). G 196-3 B 
1998ll. P Z Tel B dBiller et aklboiol ). 2M 1207 A (|Mohantv et ahlboOTl v and TWA 5B 



19991 ). In addition, we plot JHK photometry for a list of spectroscop ically con- 



firmed Taurus mem bers that do not possess circumstellar disks (iLuhman et al.l l2010l ) or known 
binary companions (jKraus et al.l l201ll . l2012l ). The JHK photometry (taken from 2MASS) has 
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been dereddened using the extinction values from the com pilation oflKraus &: Hillenbrand! ((20091) 
or the discovery survey, using the relative extincti ons from ISchlegel et al.l (|l998l ). and converted to 
MKO using the conversion of lLeggett et al.l (|2006l ). 



Rayner et alJ (120091 ) note that this figure bifurcates for field dwarfs and giants with M spectral 



types; and indeed, while some giants are on the lower sequence, no field M dwarfs are found on the 
upper sequence. HD 1160 B and C fall clearly on the upper branch with the M giants, and not the 
lower branch with the M dwarfs. HD 1160 B and C cannot be giant stars themselves, as this would 
be inconsistent with their absolute magnitudes. A small fraction of the Taurus objects lie on the 
upper sequence, but most lie on the lower branch. That so few of these young (~2-5 Myr) objects 
fall on the upper branch appears to rule out the possibility that young objects (with intermediate 
surface gravity between giants and dwarfs) populate the upper branch and older objects the lower 
one. Other plausible explanations for the five Taurus objects with giant-like colors include the 
presence of an undetected binary or disk, or an incorrect reddening determination. Or, perhaps 
some small fraction of young stars lie on the upper branch, while most lie on the lower branch. 
Alternatively, the discrepancy could also be accounted for by an error in our photometric reduction: 
shifting our NICI photometry of HD 1160 B and C by ~2cr would also place both objects on the 
lower track. In any case, with our current data we are unable to explain why the JHK photometry 
of HD 1160 B and C diverges from the photometry of dwarfs and the bulk of the single-star diskless 
Taurus objects. 

To assign a spectral type for H P 1160 C, we plot our IRTF/SpeX spectrum against M-dwarf 
standards from lCushing et ah ( 2005 ) in FigurelH We match flux levels for the entire JHK spectrum 
using only the H and K portions of the spectra, since our J-band data is particularly noisy and 
likely contaminated by HD 1160 A. The best match to the shape of the continuum and spectral 
features comes from the M3 and M4 templates. As a result, we assign a spectral type of M3.5 it 
0.5 to HD 1160 C. 



We have also compared the spectrum of HD 1160 C to an M3.5 III spectrum from lCushing et al.l 



(|2005l ) ■ but find a poor fit to the IRTF/SpeX spectrum, for both the shape of the continuum and 
strength of the absorption features. This is as we would expect: HD 1160 C, while having NIR 
colors similar to a giant, is not a giant itself. 

Figur e 1101 compares the e quivalent widths of atomic features for HD 1160 C w ith those for 



stars from lRavner et al.l (I200E) and members of the ~10 Myr T W Hydra Asso ci ation (j Covey et al 



2010l ) . As in iRayner et al.l (l2009l ) , we follow the method of ICushing et al.l (l2005l ) to compute 
equivalent widths and associated errors (see their Section 3.4), by defining continuum regions, 
fitting a polynomial, and directly computing the flux difference between measured spectrum and 
estimated continuum in the line regions. A spectral type of M3.5V is consistent with HD 1160 C. 
We cannot use these data to constrain the age, since at a spectral type of M3.5, there is not a 
significant offset between the field M dwarfs and young M stars in TWA in the strength of the 
absorption lines. 
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3.3. Age of the HD 1160 System 



Determining the age of the HD 1160 system is challenging, since the most widely used age 
diagnostics (such as lithium absorption, calcium emission. Ha emission, and ro tation rate) are 
suitable only for stars with spectral types of F, G, and K (e.g. ISoderblomI l20ld ) . In the case of 
HD 1160 ABC, with spectral types A, L, and M, none of the components are amenable to these 
techniques for determining ages, and so we turn to an analysis of the system's position on the HR 
diagram to estimate its age. In fact, this is a more fundamental method of determining ages, since 
lithium and other indirect age indicators are calibrated by measurements from open clusters, whose 
ages are determined by comparison to isochrones on the HR diagram. 

HD 1160 A was originally selected as a target for the Gemini NIC I Plan et-Fi nding Campaign 



due to its faint absolute magnitude for its spectral type. lJura et al.l (jl998l ) and iLowrance et al 



(j200d ) note that on the HR diagram early-type stars in young clusters t end to be lower (i.e. have 
smaller luminosities) than stars in older clusters. In fact. ISu et al.l (120061 ) flag HD 1160 A as lying 
below the zero-age main sequence, implying an age of <5 Myr. ( HR 8799, the host of 4 directly 
imaged planets, was determined to be young by the same method [Marois et al.ll2008l |.) With our 
discovery of two co- moving companions to HD 1160 A, we can expand this HR diagram analysis 
to include all three components, with spectral types ranging from AO to early L. 



^Figure [TT] shows HD 1160 A and C compared to stars at six sets of ages: Upper Sco (IKraus &: Hillenbrand 

20081 ) at 5 Myr; members of ~10 Myr (J3 Pic and TW Hya) and ~30 Myr (Tuc/Ho r, Carina, and 



Columba) moving groups (ITorres et al.l 12008 ) ; the 120 Myr Pleiades open cluster (IStauffer et al 



20071 ): the Hyades open cluster (iRoser et al.ll201ll ) at 600 Myr; and low- mass field objects from 
the Palomar/MSU Nearby-Star Spectroscopic Survey (iHawley et al.lll996l ). All NIR photometry 



of reference objects were obtained with 2MASS filters, and so we use the conversion factors of 
Leggett et al.l pOod ) to convert the HD 1160 ABC MKO magnitudes into 2MASS magnitudes for 
Figure [TT] (we do the same for Figure [T2]) . As we currently lack spectra for HD 1160 B and cannot 
determine its spectral type, we do not plot it in Figure [TTJ From these plots, the HD 1160 system 
appears younger than the Hyades (given the low luminosity of HD 1160 A with respect to the 
Hyades A star sequence), but consistent with the younger associations. 

In Figure [T2\ we show a NIR color-magnitude diagram for the same associations, along with 
all three members of the HD 1160 system. Here, HD 1160 again appears younger than the Hyades 
and the field objects, and inconsistent with the Pleiades as well, with HD 1160 A below the A-star 
sequence, and HD 1160 B and C redder than M stars of similar absolute J magnitudes. All three 
components are consistent with similar objects in Upper Sco or the 10 and 30 Myr moving groups 
(though the paucity of known low-mass members to the two sets of moving groups makes a detailed 
comparison at these ages difficult). As a result, we set an upper limit for the age of the HD 1160 
system of 100 Myr, and adopt an age range of 50]']4q Myr. 
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3.4. Metallicity 



An alternate explanation for the underluminosity of HD 1160 A is that the system has asub- 
solar metallicity, as we illustrate in Figure [T3l While HD 1160 A lies below the lSiess et alj (|200d ) 
tracks at [Fe/H] = 0.0, its position on the color-magnitude diagram is co nsistent for ages up to 300 
Myr at [Fe/H] = —0.3, half the metal abundance of the Sun. Similarly, iMerm et al.l (120041 ) found 
that the low position of HD 141569 (a pre-main sequence B9.5V with a circumstellar disk) on the 
HR diagram as well as its medium-resolution optical spectrum could be explained by a metallicity 
of [Fe/H] = -0.5. 

We note, however, that Figure [T2] provides two indications of youth for HD 1160: the under- 
luminosity of A and the redness of B and C in J — Ks color. In Figure Owe sho w a NIR color- 



magnitude diagram with HD 1160 B and C and field M dwarfs of known metallicity (jLeggett et al 



2000 



Roias-Avala et al.ll2010l ). M dwarfs with subsolar metallicity have bluer J — Kg colors than 



more metal-rich M dwarfs. So if the HD 1160 system were to be old (>100 Myr) but metal-poor, 
we would expect bluer J — Ks colors for HD 1160 B and C than we observe. The underluminosity 
of A stars on color-magnitude diagrams may indicate either young age or low met allicity, but the 
only a vailable explanation for red J — Kg colors in late-type stars is youth. Indeed, iJohnson et al 



(j201ll ) derive a relationship between J — Kg color and metallicity of M dwarfs, with redder J 
color indicating more metal-rich atmospheres. 



Ks 



We also note that regardless of whether HD 116 A has [Fe/H] = 0.0 or —0.3, the expected 
main sequence lifetime of an AO star is ~300 Myr (ISiess et al.l 12000) , which provides an upper 
limit on the age of the system independent of the measured photometry of HD 1160 B and C. 
Additionally, since st ars recently formed in the solar neighborhood tend to be of solar metallicity 
([Pedicelli et al.ll2009l ). for an age less than ~300 Myr, the HD 1160 system would have had to 
have formed further out in the galactic disk (i.e. in a lower metallicity environment) and migrated 
inward, which is inconsistent with its main sequence lifetime of ~300 Myr and its low space motion 
(see Section [3.5p . Though a metallicity of [ Fe/H] = —0.3 is at the lower limit of Cepheids at the 
solar radius, as measured bv lPedicelli et al.l ( 20091 ). the more likely explanation is that HD 1160 A 
is solar metallicity, and its low position in the HR diagram is due to youth. 



3.5. Space Motion 



By combining our measured RV of the HD 1160 system w ith the position and revised Hipparcos 
parallax and proper motion of HD 1160 A (Ivan Leeuwenll2007l ). we have computed the space motion 
of HD 1160. We find the space motion to be C/ = -7.6 ± 0.4, V = -3.4 ± 0.5, W = -15.7 ± 
0.4 km/s. In Figure 1151 we co mpare this space motion to those of young, nearby moving groups, 
as given bv lTorres et al.l (j2008l ). HD 1160 is not co-moving with any known moving group, though 
its space motion is relatively low, with a total motion of 17.7 it 0.5 km/s. Such a low velocity is 
consistent with a young age for HD 1160. We also note that this UVW motion of HD 1160 places 



-li- 



the system s lightly outsid e (though near the edge) of the la UVW ellipse of the "young disk" 
computed bv lEggenl (jlOsd l oi U = -15 ± U, V = -14 ± 9, = -6 ± 13 km/s. 



3.6. Mass of HD 1160 B and C 



Using the age and available photometry and spectroscopy of HD 1160 B and C, we pro- 
ceed to assign masses to these companions using theoretical evolutionary models for stellar and 
sub-stellar objects. Table H] s hows mass estiinates from the Lyon evo lutionary models using the 
DUSTY atmospheric models (jchabrier et ahlboool : lAllard et ahlboOlh for HD 1160 B for the ab- 
solute magnitudes in each infrared band. Errors are calculated by interpolation across gridpoints 
of the DUSTY models, and using a Monte Carlo analysis assuming gaussian distributions for the 
magnitude uncertainties and a uniform distribution of ages from 10-100 Myr. Given the large age 
range, the implied mass range is quite large (24-90 Mj^p), spanning much of the range of brown 
dwarf masses. We note that these ranges only account for uncertainties in fluxes and the age, and 
do not account for any systematic uncertainty in the models themselves. From the J-band (the 
peak of the NIR SED for late M and early L spectral types) flu x alone, w e derive a mass for HD 
1160 B of SSlg^ Mjup. Using the BC( J) bolometric correction of IlIu et al.l (|2010al ) (BC( J) chosen 
since it is flat over the plausible range of spectral type for HD 1160 B, L0ib2), we find a similar 
mass by computing the bolometric luminosity (L = 5.9x10"^ it 5xl0~^ L©); cind comparing to 

+12 



the DUSTY models: M = 371^5 Mjup. 

For HD 1160 C, we incorporate our measured s pectral type of M3.5 ± 0.5 when estimating the 
mass. Using the dwarf sequence of iLuhmanI (|l999l ). we convert this spectral type to an effective 
temperature of 3270 ± 90 K, (or 3340 ± 70 K usiri g the intermediate scale) and using a K-band 
bolometric correction from iGolimowski et al.l (|2004l ). we comp ute a luminosi t y of H D 1160 C of 
9.4x10"^ ± 8x10""^ Lq. Coupled with the NextGen models of baraffe et al.l (Il998ll . these values 
allow us to estimate a mass for HD 1160 C of I90I40 ^Jup- Using the |Luhman| |l99i) intermediate 



sequence (more appropriate for younger stars), we find the mass of HD 1160 C to be 230^45 ^Jup 
(0.22l0;03 Me). 



3.7. Orbital Stability 

We estimate the likely sen iimaior axes o f the tw o companions from their projected separations 
using the conversion factor of iDupuy fc Liul (1201 ll ): assuming a uniform eccentricity distribution 
between and unity and no discovery bias, their estimated ratio of semi-major axis (a) to projected 
separation (p) is a/p = 1.10 with 68.3% confidence limits between 0.75 and 2.02. Given projected 
separations of 81 it 5 AU and 533 it 25 AU for HD 1160 B and C, and using a Monte Carlo method 
to propagate errors (using a gaussian distribution for separation and two gaussians to account for 
the assymetric error bars on a/p), we estimate semi-major axes of 74 it 35 AU for HD 1160 B and 
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500 ± 200 for HD 1160 C. 

The HD 1160 ABC system appears to be orbitally stable, assuming that the eccentricity of 
the A and C components is not larger than about 0.3. Given the mas s ratio between HD 1160 A 



and C of 0.1, and using the the detailed orbital stability calculations of lHolman &: WiegertI (|l999l ). 
and assuming an orbital eccentricity of these two major components of 0.3, we estimate that for 
the HD 1160 system the critical semimajor axis for stability of component B is about 0.2 times the 
semimajor axis of components A and C. That is, component B can orbit no farther from A than 
about 100 AU. With an estimated semi-major axis of about 74 AU, component B therefore appears 
to be part of a stable hierarchical system. Smaller orbital semimajor axes for HD 1160 B would 
also be stable, and so additional companions could plausibly exist interior to B. 

Since we have not fit an orbital solution to our astrometry, the semi-major axes of both 
components may differ from our estimates. Indeed, if the semi-major axis is 109 AU for HD 1160 B 
and 300 AU for HD 1160 C (the la upper and lower limits of their respective semi-major axis 
estimates), then the system would not be orbitally stable. Additionally, a high eccentricity for one 
or both of HD 1160 B and C would suggest instability. Nevertheless, given our measurements, and 
the small probability of catching this system in an unstable state, the most likely explanation is 
that the system is stable. 



Conclusions 



We have described here the discovery and analysis of two low-mass companions to the young 
A star HD 1160 A, found during the Gemini NICI Planet-Finding Campaign. With estimated 
companion masses of 33lg^ and 230^*^45 ^ J^p ^^'^ ^ primary mass of ~2.2 Mq, this system has very 
small mass ratios between components, 0.10 between A and C, and 0.014 between A and B. Around 
a star slightly less massive than the Sun, this A-to-B mass ratio would represent a planetary-mass 
companion. 

Hierar chical triple systems li ke HD 1160 occur with a frequency of about 8% around solar- 
type stars (jRaghavan et al.ll2010l ). though seldom with such extreme mass ratio s as the A to B 



compc 
20l3) 



components of HD 1160. Mass ratios less than about 0.1 are rare (see Figure 16 of lRaghavan et al 
tho ugh not unpreceden ted, e.g., ( Vir B, an M-dwarf companion to an A star, with mass 



ratio 0.08 (jHinklev et al.ll20ld ) resembles HD 1160 C, while the brow n dwarf HR 7329 B, with mass 
ratio of ~0.01 relative to its A-star primary (jLowrance et al.ll200d ). is more akin to HD 1160 B. 
Along with binary star systems, multiple star systems like HD 1160 are thought to form by the 
collapse and fragmentation of dense molecular cloud cores (e.g., lBossll2009l and references therein). 
The cloud core that collapsed to form HD 1160 presumably was relatively rapidly rotating, in order 
to result in the wide orbits of all three components. Typical separations for binary stars with mass 
ratios of about 0.1 ( as for A and C) are in the range of about 10 to 10^ AU (see Fi gure 17 of 
Raghavan et al.ll2010l ). so a separation of A and C of about 531 AU is consistent with this range. 
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Near-IR spectra of HD 1160 B will likely set much stronger constraints on the age of the 
system, as an ~L0 brown dwarf (as suggested by the photometry shown in Figure [7|) will show 
significant spectral evolution up to the main sequence lifetime of HD 1160 A (^300 Myr). If the 



spectrum shows s i gns o f low surface gravity compared to field L dwarfs (e.g. iMcGovern et al.ll2004l : 



Allers et al.l 120071 . l2009l ). this will help refine our age determination of the system. 



In the longer term, high-precision orbital monitoring of both components may be able to 
constrain the orbital parameters of HD 1160 B and C. Given the expected orbital periods of ~400 
and ~7000 years, such orbital motion will be quite small and so will require high precision and 
long time baselines. In fact, the orbital period of HD 1160 B may be significantly longer than 400 
years, since a circular face-on orbit would predict ~120 mas of motion over the 10 years for which 
we have archival data, and a linear fit to the relative motion of HD 1160 B gives only 23 ± 29 mas 
over 10 years (< la, so not a significant detection of orbital motion). The fact that we observe a 
larger degree of orbital motion from the further-out C (60 zt 30 mas over 10 years) than from the 
close-in B suggests that there may be significant misalignment of the orbital planes of these two 
companions. Additional follow-up and analysis of this system will set interesting constraints for 
modeling the formation of HD 1160 ABC. 
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Fig. 1.— Images of HD 1160 A, B, and C over twelve epochs from 2002 to late 2011. The VLT 
NACO data were taken as calibration data for science programs targeting other objects, and in 
some cases HD 1160 B is detected at low S/N. Each image was rotated to place North up and East 
to the left and has the same field of view, 2.7" x 6.5". 
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Fig. 2.— Archival images of HD 1160 A and C from 2002 until 2011, taken with the VLT ISAAC 
instrument. As with the VLT NACO data, these are calibration data and are sometimes at low 
S/N. Additionally, since the seeing- limited resolution of ISAAC is lower than the AO-corrected 
images in Figure [H HD 1160 B is not detected at any of these epochs. Each image was rotated to 
place North up and East to the left and has the same field of view, 2.7" x 6.5". 
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Fig. 3. — Astrometry of HD 1160 B from 2002 to 2011 (points with error bars) along with the 
expected motion of a background object (black line). A background object should fall within the 
gray shaded track, given the known proper motion and parallax of HD 1160 A. The thickness of 
the gray track indicates the positional uncertainties arising from the proper motion and parallax 
uncertainties for HD 1160 A, as well as the astrometric uncertainties from the NICI 2010.8301 
position to which the track is tied. HD 1160 B is clearly a common proper motion companion. 
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Fig. 4. — Astrometry of HD 1160 C from 2002 to 2011. While a less significant detection of common 
proper motion than with HD 1160 B, the data are still much more consistent with HD 1160 C being 
a co-moving companion rather than a distant background object. 
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Fig. 5. — Astrometry for HD 1160 B with respect to HD 1160 A as seen on the sky. Each data 
point is connected by a thin hne to the expected location of a distant background object at the 
observational epoch. The labels on the right give the epoch of the astrometric measurements, with 
the vertical location of the epoch label set by the Declination expected for a distant background 
object at the epoch. Within the error bars, all the datapoints are consistent with the position mea- 
sured at the NICI 2010.8301 detection and are significantly displaced from the expected positions 
for a distant background object. The reduced value for HD 1160 B being a background object 
is 26.2 (dof=22, Pf«0%), compared to xt = 0-38 (dof=22, P=99.6%) for a common proper motion 
companion. 
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Fig. 6. — On-sky astrometric plot for HD 1160 C with respect to HD 1160 A. The fit to the 
motion of a common proper motion companion is much more hkely than the fit to the motion of 
a background object: the significance of the fit is Xu — ^-9^ (dof=36 P~ 0%) for a background 
object and xt ~ 0-62 (dof = 36, P=96.5%) for a common proper motion companion. In addition, 
the expected direction of the motion of a background object is orthogonal to the spread of the 
datapoints, and we have measured the radial velocities of HD 1160 A and C to be consistent within 
measurement uncertainty. 
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Fig. 7. — Color-color diagram on the MKO system, suggesting spectral types of early-L and mid-M 
for HD 1 160 B and C, respec tively. The low-mass field dwarfs plotted here are from the compi- 
lation of iLeggett et al.l ([20101). The implied spectral type for HD 1160 C is a close match to our 
determination of M3.5 based on near-IR spectroscopy. 
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Fig. 8. — Color-color diagram showing HD 1160 B and C, along with giants, dwarfs, single-star 
members of Tam'us, and some individual low-mass objects. While HD 1160 B and C are not giants 
themselves, they do exhibit giant-like JHK colors. 
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Fig. 9. — IRTF/Sp eX spectra of HP 11 60 C compared to M dwarf spectral standards and line 
identifications from ICushing et al.l (j2005l ). The best fits to the shape of the H and K continuum 
are spectral types M3 and M4, and therefore we assign a spectral type of M3.5. Note that the 
J-band portion of the spectrum suffers from contamination from HD 1160 A (AO), and so we only 
use the H and K portions in our fits. 
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Fig. 10. — Comparison of equiv alent widths o f atomic features in the HD 1160 C spectrum to those 
of dwarf and giant stars fr om iRavner et al.l (|2009l ) and members of the TW Hydra Association 
(TWA) from ICovev et al.l (j20ld ). Red squares denote dwarfs, green triangles giants, and black 
crosses TWA members. For each feature, the black horizontal line denotes the equivalent width of 
HD 1160 C and the gray swath the corresponding measurement uncertainty. Taken together, the 
atomic absorption lines support our M3.5 spectral type for HD 1160 C. However, these data cannot 
constrain the age of HD 1160 C, since the equivalent widths for M dwarfs and TWA members are 
similar. 
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Fig. 11. — HR diagram of stars in Upper Sco, 10 Myr moving groups (/3 Pic and TW Hya), 30 
Myr moving groups (Tuc/Hor, Carina, and Columba), the Pleiades, the Hyades, and field low-mass 
objects (magnitudes are in the 2MASS system). The positions of HD 1160 A and C are marked 
with filled red circles. We do not include HD 1160 B, as we do not have a spectral type. HD 1160 A 
appears underluminous compared to other A stars, especially for older stars in the Hyades. 
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Fig. 12. — Color-magnitude diagram (2MASS system) showing HD 1160 ABC with respect to stars 
in stellar associations of different ages. HD 1160 A is too underluminous to fit either the Hyades 
or Pleiades stars, while HD 1160 B and C are too red to fit the low-mass end of the main sequence 
of these older open clusters or field objects. The HD 1160 system does appear consistent with the 
stars in Upper Sco and the 10 and 30 Myr moving groups. As a result, we assign an age to the 
system of bOtlo ^y^- 
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Fig. 13. — Optica l color-magnitude diagram showing the position of HD 1160 A (filled circle) against 



Siess et al.l ((20001) isochrones, for solar metallicity (solid lines) and half-solar abundance (dashed). 
HD 1160 A is underluminous compared to tracks of all ages at [Fe/H] = 0.0, but is consistent with 
ages of 10-300 Myr for [Fe/H] = —0.3. However, the colors of the B and C components do not 
agree with a low metallicity for the system (see Figure [1 
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Fig. 14. — HD 1 160 B an d C compared to M dwarfs of known metallicity, compiled by 



Roias-Avala et al.l (|2010l ) and iLeggett et al.l (|200d ) . Low-metallicity M dwarfs tend to be bluer 
in J — i^T color than solar-metallicity dwarfs, especially at later spectral types. So while the under- 
luminosity of A stars such as HD 1160 A can be explained by either youth or sub-solar metallicity, 
the J — K redness of HD 1160 B and C supports the idea that the system is young. 
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Fig. 15. — UVW motions of nearby young, moving groups from [Torres et al.l (j2008l ) compared to 
the HD 1160 system (BP: /3 Pic, TH: Tuc/Hor, Col: Columba, Car: Carina, TWA: TW Hya, eC: 
rj Cha, Oct: Octans, Arg: Argus, ABD: AB Dor). Error bars for the moving groups indicate the 
RMS of the space motions for the association members, while the error bars for HD 1160 indicate 
measurement uncertainties. HD 1160 is not co-moving with any well-known association, though 
the magnitude of its space motion is consistent with these young stars, supporting a young age 
(<100 Myr) for HD 1160. 
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Table 1. Properties of HD 1160 A 



Property 



Measurement 



Spectral Type 
RA (ep J2000) 
Dec (ep J2000) 
Proper Motion RA (mas/yr) 
Proper Motion Dec (mas/yr) 
Parallax (mas) 
Radial Velocity (km/s) 
U (km/s) 
V (km/s) 
W (km/s) 
J (mag) 
H (mag) 
Ks (mag) 
L' (mag) 
Ms (mag) 
Age (Myr) 



AOV 
00:15:57.3025 
+04:15:04.018 
21.15 ± 0.62 ^ 
-14.20 ± 0.24 ^ 
9.66 ± 0.45 ^ 
12.6 ± 0.3 ^ 
-7.6 ± 0.4 ^ 
-3.4 ± 0.5 ^ 
-15.7 ± 0.4 
6.98 ± 0.02 
7.01 ± 0.02 = 
7.04 ± 0.03 = 
7.055 ± 0.014 
7.04 ± 0.02 <^ 



50 



+50 b 
-40 



van LeeuwenI (|2007l ) 



'^This work 



Cutri et al.1 (120031 1 



' iLeggett et all ([2003) 
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Table 2. Astrometry of HD 1160 BC 



HD 1160 B HD 1160 C 







Sep ( 


J 




PA (°\ 




Sep ( 


J 




PAr; 


) 




lllbLl Ulllclll 


2002.5699 


0, 


.77 




0. 


,03 


246.2 ± 1.0 


5, 


,17 




0. 


,03 


350.9 




0. 


,5 


VLT NACO Ks 


2002.7397 














5, 


.17 




0, 


,04 


349.8 




0, 


,5 


VLT ISAAC L 


2002.8658 














5, 


.22 




0, 


,04 


349.6 




0, 


,5 


VLT ISAAC L 


2003.5342 














5, 


.19 




0, 


,04 


349.5 




0, 


,5 


VLT ISAAC L 


2003.6356 














5, 


.18 




0, 


,04 


349.7 




0, 


,5 


VLT ISAAC L 


2003.8384 


0, 


.77 




0, 


,03 


245.6 ± 1.0 


5, 


.16 




0, 


,03 


350.1 




0, 


,5 


VLT NACO Ks 


2003.8603 














5, 


.18 




0. 


,04 


349.6 




0. 


,5 


VLT ISAAC L 


2005.9753 


0, 


.76 




0. 


,03 


244.7 ± 1.0 


5, 


,15 




0. 


,03 


350.4 




0. 


,5 


VLT NACO V 


2007.9014 














5, 


.08 




0. 


,07 


349.4 




0. 


,5 


VLT ISAAC K, 


2008.5027 


0, 


.80 




0. 


,06 


245.3 ± 2 


5, 


,16 




0. 


,03 


349.7 




0. 


,5 


VLT NACO L' 


2010.7096 


0, 


.77 




0. 


,06 


242.8 ± 2 


5, 


.15 




0. 


,03 


349.4 




0. 


,5 


VLT NACO V 


2010.8301 


0, 


.78 




0. 


,03 


244.3 ± 0.2 


5, 


.15 




0. 


,03 


349.8 




0. 


,2 


Gemini NICI 


2010.8904 


0, 


.76 




0. 


,03 


244.5 ± 0.2 


5, 


.16 




0. 


,03 


349.6 




0. 


,2 


Gemini NICI H 


2010.9041 


0, 


.78 




0. 


,02 


244.9 ± 0.5 


5, 


,14 




0. 


,02 


349.9 




0. 


,5 


Keck NIRC2 V 


2011.5233 


0, 


.78 




0, 


,03 


244.0 ± 1.0 


5, 


.14 




0, 


,03 


349.4 




0, 


,5 


VLT NACO L' 


2011.6685 


0, 


.78 




0. 


,03 


244.9 ± 1.0 


5, 


.14 




0. 


,03 


349.5 




0. 


,5 


VLT NACO V 


2011.7589 














5, 


.12 




0. 


,07 


349.2 




0. 


,5 


VLT ISAAC /C, 


2011.8027 


0, 


.77 




0. 


,03 


244.5 ± 0.2 


5, 


.16 




0. 


,03 


349.6 




0. 


,2 


Gemini NICI H 


2011.8521 


0, 


.78 




0. 


,03 


244.4 ± 1.0 


5, 


.14 




0. 


,03 


349.4 




0. 


,5 


VLT NACO V 



-37- 



Table 3. Photometry of HD 1160 BC 



Bandpass 


HD 1160 B 


HD 1160 C 


Instrument 


(MKO) 


(mag) 


(mag) 




AJ 


8.85 ± 0.10 


6.33 ± 0.04 


Gemini NICI 


AH 


7.64 ± 0.08 


5.53 ± 0.03 


Gemini NICI 


AKs 


7.08 ± 0.05 


5.14 ± 0.06 


Gemini NICI 


AL' 


6.35 ± 0.12 


4.803 ± 0.005 


Keck NIRC2 


AMs 


7.3 ± 0.2 


5.10 ± 0.05 


Keck NIRC2 


Mj 


10.75 ± 0.14 


8.23 ± 0.11 




Mh 


9.57 ± 0.13 


7.46 ± 0.11 




Mks 


9.04 ± 0.12 


7.10 ± 0.12 




Ml, 


8.33 ± 0.16 


6.78 ± 0.10 




Mms 


9.3 ± 0.2 


7.06 ± 0.11 
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Table 4. Mass Estimates for HD 1160 BC 



Parameter Measurement Mass 

(Mj„p) 



HD 1160 B 



Mj (mag) 


10.75 ± 0.14 


33+^^ 


a 


Mh (mag) 


9.57 ± 0.13 


48lJ^ 


a 


Mks (mag) 


9.04 ± 0.12 




a 


Ml' (mag) 


8.3 ± 0.16 


59^?J 


a 


Mms (mag) 


9.3 ± 0.2 


qq+20 


a 


BCj 


2.06 ± 0.14 ^ 


07+12 

-12 


a 



HD 1160 C 



Sp. Type M3.5 ± 0.5 
BCa- 2.72 ± 0.06 

Dwarf Temperature Scale 
Te// 3270 ± 90 1901 

Intermediate Temperature Scale 
Te// 3340 ±70'^ 2301 



Chabrier et al.1 tood ) 



Uiu et al.1 (|2010al l 



Golimowski et al.1 (|2004l ) 



^Luhmanl (ll999|) 



Baraffe et al.l (|l998l ^ 



